Does the enhancement observed in yy — » DD contain two P-wave higher charmonia? 
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Solved is a new puzzle raised by the observation of an enhancement structure Z(3930) in yy — > DD. If cate- 
gorizing Z(3930) slsXci(2P) suggested by Belle and BaBar, we must explain why Xco(2P) dominantly decaying 
into DD is missing in the DD invariant mass spectrum. In this letter, we propose that the Z(3930) enhancement 
structure may contain two P-wave higher charmonia XcO (2/*) and / ^ c j(2P). We show that this assumption is sup- 
ported by our analysis of the DD invariant mass spectrum and cos 6* distribution of yy — > DD. In addition, we 
also extract the resonance parameters of ^ d) (2P) and Xdi^P) which show that these are two narrow resonances 
with masses very close to each other, i.e., M XcB (2P) = 3.92 GeV with T xm q P) = 8.07 MeV and M Xc1 q P) = 3.94 
GeV with T Xc2i2P> = 11.98 MeV. These results would not only provide valuable information of two P-wave 
higher charmonia Xco(2P) and;^ C 2(2P), but also serve as the crucial test of our novel proposal to the observed 
enhancement structure Z(3930), especially at the forthcoming Bellell and the approved SupreB. 

PACS numbers: 13.25.Gv, 12.38.Lg 



Introduction : Since the Belle Collaboration reported the 
first observation of X(3872) [1], more and more charmonium- 
like states XYZ have been observed in experiments [2]. 
Both theoretical and experimental studies of charmonium- 
like states have become an intriguing and challenging re- 
search field of hadron physics. It is well known that the 
potential model of hadron has made success for describing 
the low-lying hadron spectra [3] as well as heavy quarko- 
nia. Meanwile, the abundant higher charmonium and bot- 
tomonium spectra were predicted. We believe these newly ob- 
served charmonium-like states have a close relationship with 
higher charmonia, which provides a good platform to per- 
form the study of higher charmonium. It is valuable to reveal 
not only the underlying structures of these charmonium-like 
states, but also our deep understanding of non-perturbative be- 
havior of QCD. 

As listed in Particle data group (PDG) [4], there are 
three established P-wave spin-triplet charmonia, ^ f o(3415), 
^f c i(3510), and ^ C 2(3556). However, compared with the S- 
wave, the P-wave charmonium spectrum has not yet been 
established. The observations of charmonium-like states 
X(3812) [1], Z(3930) [5, 6] and X(3915) [7] also have stim- 
ulated the theorists' extensive interest [8-10] in exploring the 
higher excitations of these P-wave spin-triplet charmonia. 

Together with these observations, we need to know details 
of Z(3930) [5, 6], which is a charmonium-like state observed 
by Belle [5] and confirmed by BaBar [6] in the yy — > DD pro- 
cess. Belle gave its mass M = 3929 ± 5(stat) + 2(syst) MeV 
with T = 29 + 10(stat) ± 2(syst) MeV [5], and BaBar obtained 
the consistent results with Belle, M = 3926.7 + 2.7(stat) ± 
l.l(syst) MeV with T = 21.3 + 6.8(stat) ± 3.6(syst) MeV [6]. 
Both Belle and BaBar indicated that Z(3930) is the candidate 
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of Xcii^P) suggested by their analysis of the angular distri- 
bution [5, 6]. Thus, the Z(3930) experimentally established 
makes the P-wave spin-triplet charmonium spectrum abun- 
dant. In addition, it is not difficult to explain X(3872) as the 
first radial excitation of Xc\ (35 10) when considering the cou- 
pled channel effects [8, 10] or the dominant 2P cc component 
mixing with D°D*° + D M D° [11]. 

As the P-wave spin-triplet charmonium spectrum becomes 
more abundant, however, an urgent and crucial question 
emerges out of the study on the first radial excitation of P- 
wave charmonia. The mass of ^- c o(2P) theoretically predicted 
to be the first radial excitation of ^ f o(3415) is very close to 
that of Z(3930) [8, 10] and above the DD threshold. Addi- 
tionally, Z(3930) decays into DD via the D-wave interaction, 
while ^- c( )(2P) via the S -wave, where ^f e o(2P) — > DD is a dom- 
inant contribution to the total width [9]. Since Z(3930) was 
already observed in the DD invariant mass spectrum of the 
yy — > DD process [5, 6] and we cannot find any suppression 
mechanism of the ^ c o(2P) production in this process, where 
^ c o(2P) and Z(3930) are assumed to have the same spatial 
wave functions, we believe that ^ c o(2P) should exist in the 
same data samples of the DD invariant mass spectrum. How- 
ever, the present experiment did not report any evidence of 
^ e o(2P) in this process [5, 6], which obviously contradicts the 
above general analysis. This is a new puzzle when studying 
the P-wave higher charmonia. 

Due to the peculiarity of x c q(2P), in this letter we propose 
a novel conjecture to solve this puzzle, i.e., the enhancement 
structure observed in the DD invariant mass spectrum of the 
yy — > DD process [5, 6] should contain both ^ f o(2P) and 
XciQP) signals. To testify that this conjecture is reasonable, 
in the following analysis we shall construct a model depicting 
the yy — > DD process, where we will consider both^ c o(2P) 
and XciQP) contributions to this process. Furthermore, by 
fitting the experimental data of the DD invariant mass spec- 
trum [5, 6], we will test this conjecture and extract the reso- 
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nance parameters of Xco(2P) and;^2(2P), which are important 
properties of P-wave higher charmonia. What is more impor- 
tant is that we find a complete series of P-wave spin-triplet 
charmonia including the established ground states and their 
corresponding first radial excitations. Of course, this will in- 
spire experimentalists' interest in studying these Xco(2P) and 
Xa{2P) in future experiment, too. 

Formulation : In general, 7(^1)7(^2) — > D(p\)D(p2) oc- 
curs via two different mechanisms. The first mechanism is the 
so-called direct process, which provides the background con- 
tribution when studying the DD invariant mass spectrum. Its 
amplitude can be parameterized as 



where 




D 

1 



D 



8n. 



oR ^e v 2 (T^ + c 02 T^)r(s),(l) 



which is abbreviated as ZRnor- 
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£>fiv,aj3 , P a ^ pA pipP2A correspond to the helicity-0 and helicity-2 
contributions, respectively, which can be distinguished by the 
superscripts and 2. Here e^(ki) and e^iki) are polarization 
vectors of two photons and the definitions of 7^ v , & MV ,a/3, and 
<p a 0P A w in be given later. T{s) is the function of the energy in 
the 77 center-of-mass frame, and is defined as [6] 



-Js - m 2 { - m,Texp[-/3( - m,)] /mf + \ (2) 



where m, is taken as the DD threshold. Free parameters a 
and j6 are determined by fitting the experimental data. In or- 
der to make the form factor dimensionless, a factor mj <a+1> is 
multiplied. 

The second mechanism is the intermediate resonance state 
contribution to the 77 — > DD process, which is from the s- 
channel. The initial 77 and final DD are connected by the 
intermediate resonances with J p = + ,2 + , where the parity 
and total quantum number of the intermediate state are con- 
strained by both the initial 77 and final DD. When analyzing 
the data of Belle [5] and BaBar [6], we select Xco(2P) and 
Xci(2P) as the intermediate states, which are the first radial 
excitations of ^^(3415) and ^ f 2(3556), respectively. These 
ground states can also contribute to yy — » DD as the interme- 
diate resonances. However, since the masses of^- ri )(3415) and 
^ C 2(3556) are far away from the energy range discussed here, 
the high-mass tail of these resonances only provide one of the 
backgrounds to the DD invariant mass spectrum in the energy 
range (> 3.7 GeV), which can be included in Eq. (1). 

To investigate the process yy — > DD by the intermediate 
Xm{2P) and Xci{2P), we adopt the effective Lagrangian ap- 
proach. Here, the interactions of 77 and DD with Xcj(2P) are 
described by the effective Lagrangians 

£ X J>(2P)yy ~ gxco(2P)yyA , 'A y T fl vXco(2P), (3) 

■QaOPyrr = gx,7(2P)y y A' 1 A v x a C 2{2P)& llv , a p, (4) 

-C Xc0 (2P)DD = g XcB (2P)DBXco(2P)DD, (5) 

£ Xc2 (2P)dd = -g Xc2 (2P)DDX^ V 2 (2P)dfiDdvD, (6) 



T» v = g^-k v X/h-k 2 , 



[g^k^k^ + g^k^ - g^k^] 
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in the momentum space. 

The propagators of Xcj(2P) (J = 0, 2) are in the form, 



Xco(2P) : 



Xc2(2P) : 



' m XM (2P) + im xM2P)^xM2P) ' 
j tpiiviiV 



I 1 ~ m L(2P) 



withP^v = ^dw'^vv+lp/'lv/xO-j^Vv and^ v = -g^+ 

In this mechanism, one obtains the amplitudes for 77 — > 
DD wmxco(2P)andxc2(2P), 
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? ~ Xci(2P) + im Xc2(2P)T Xc2 (2P) 

The total amplitude of 77 — > DD is written as 

Mrotal = ^NOR + e®°M Xc0 (2P) + (10) 

where the background and intermediate Xc()(2P) /Xci(2P) res- 
onance contributions are considered. Phases 0o and fa are 
introduced to describe the the interference of the amplitudes. 
Finally, one can compare the measured distributions with 



— =N — =N — —\M T f 

dt dt 64ns \pi cm \ 2 



(11) 



where s = q 2 = {k\ + k-ij 2 is the square of the energy, p lcm 
is the momentum of D in the center-of-mass frame and No a 
normalization factor. In addition, we define 



■ Pif = m% - I + - yls(s-4m 2 D ) cos 0*. 



And 6* denotes the angle between outgoing D meson and in- 
coming photon in the center-of-mass frame. 

Numerical Calculation : As presented in Refs. [5, 6], 
the experimental data of the DD invariant mass spectrum and 
cos#* distribution of 77 — > DD were released. Based on 
these data, with the help of Eq. (11) we study whether this 
enhancement structure observed by Belle and BaBar contains 
both^ c ()(2P) and^ C 2(2P) resonances. 



3 



TABLE I: The values of these parameters obtained by fitting the 
Belle data [5]. The resonance parameters masses and decay widths 
°f Xd>(2P) and^ C 2(2-P) are in units of GeV and MeV, respectively. 
The phase angles, a and f Xcl (iP) are dimensionless. p and f Xda Qp) are 
in units of GeV 1 and GeV 2 , respectively. 



Parameters 


Values 


Parameters 


Values 


m Xd>(2P) 


3.920 ± 


0.007 


r Xc0(2P) 


8.065 ± 


9.663 




3.942 ± 


0.003 


r Xd(2P) 


1 1.980 1 


: 6.953 




0.035 ± 


0.023 


0o (Rad) 


2.919 ± 


0.287 


fxaQP) 


2.836 + 


1.252 


<t>2 (Rad) 


2.492 + 


0.390 


a 


-0.497 d 


: 0.788 


P 


8.996 ± 


1.172 


C 2 


39.996 ± 


11.509 









In our model, the lineshape of the DD invariant mass spec- 
trum and the cos 9* distribution are determined by the follow- 
ing 1 1 parameters: 

m Xc0 (2P), T Xca &p), m Xc2{2 p), r Xc2 ( 2 P),f Xcll (2P), f Xc2 (2P), 
a, j3, c 02 , (f>o, <f>2, 

where f Xc0 (2P) = gx c o(2P)yyg Xc0 (2P)DD/gNoR and f Xe2 (2P) = 
8xci(2P)yyS Xc2 (2P)DD/ 8nc,r and the coupling constants g NoR can 
be absorbed by the normalization factor No. By fitting the ex- 
perimental data of the DD invariant mass spectrum and cos 9* 
distribution of yy — > DD [5, 6], these parameters can be de- 
termined. Among all fitting parameters, m XcJ (2P) and F Xc] (2p) 
(J = 0, 2) are important resonance parameters reflecting the 
properties of^ c0 (2P) and^ c2 (2P). 

With the above preparation, we perform a fitting of our 
model to the Belle's data [5]. Here, the reason why we choose 
the Belle's data is that its information is more abundant than 
that from BaBar [5, 6]. For instance, Belle gave the cos#* 
distribution in the 3.91 < M(DD) < 3.95 GeV region of 
yy — » DD including only the background contribution (the 
cyan histogram in Fig. 2) [5], which enables us to constrain 
co2 in Eq. (1) appropriately. The cos 9* distribution of the 
background extracted by the Belle collaboration depends on 
cos 6*, which indicates that there is helicity-2 component in 
the background and this component will interfere with the 
resonances contributions from xaQP) state and such inter- 
ference is not included in the experimental analysis [5]. 

In Table I, we have listed the values of parameters corre- 
sponding to our best fit to the Belle's data. Fig. 1 shows the 
comparison of theoretical line shape of the DD invariant mass 
spectrum of yy — > DD with the experimental data. Although 
we fit our model to the Belle's data in the analysis of the DD 
invariant mass spectrum, we also make a comparison of our 
results with the BaBar's data. We need to normalize the m DD 
invariant mass spectrum distribution given by BaBar to the 
Belle's data by multiplying the factor, which is the ratio of 
the total events measured by Belle to those obtained by BaBar 
for the corresponding bins. As shown in Fig. 1, the obtained 
line shape of the DD invariant mass spectrum of yy — > DD 
can indeed describe these experimental data. In Fig. 1, we 
also show the line shape of the background (grey solid curve) 
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FIG. 1: (color online). The best fit (red histogram) to the experi- 
mental data of the DD invariant mass distributions given by Belle [5] 
(blue dots with error bar) and BaBar [6] (green triangles with error 
bar). Here, the grey solid curve is the description of the background 
of 77 — > DD with these obtained parameters listed in Table I and Eq. 
(1). 



when keeping only J\.nor contribution and adopting the opti- 
mum values of a, [3, and C02 listed in Table I as the input. 

We need to emphasize that a very steep line shape at 
Mod ~ 3.95 GeV exists in the experimental data of Belle 
[5] and BaBar [6], which is reflected in our fit. When in- 
cluding Xco{2P) and Xc2(2P), we can also reproduce the en- 
hancement structure in the DD invariant mass spectrum of 
yy — » DD [5, 6]. Especially, the extracted resonance parame- 
ters ofXco(2P) and^- C 2(2P) indicate that these are two narrow 
resonances with masses very close to each other. These opti- 
mum resonance parameters are also consistent with the theo- 
retical predictions of their masses [8, 10] and widths [9]. 

Figure 2 shows our fit to the cos 9* distribution of yy — > 
DD given by Belle [5]. Using our model with the parame- 
ters in Table I, we further get the cos 6* distribution in the 
3.91 < M(DD) < 3.95 GeV region of yy -> DD correspond- 
ing to Fig. 1. Comparing our results with the Belle's data, one 
finds that the data on cos 0* released by Belle can be well de- 
picted even if the enhancement structure in the DD invariant 
mass spectrum is formed by two f-wave higher charmonia 
Xco(2P) and Xc2(2P). In Refs. [5, 6], both Belle and BaBar 
once claimed that the enhancement structure observed in the 
DD invariant mass spectrum is due to only one resonance with 
J PC = 2 ++ as the candidate of ^ C 2(2P) by performing the an- 
gular distribution analysis [5, 6], while we claim that two res- 
onances can explain the same enhancement structure as shown 
by detailed numerical calculations in this letter. 

Summary : An enhancement structure named as Z(3930) 
in yy — > DD reported by Belle [5] and confirmed by BaBar 
[6] raises a new puzzle in the investigation of P-wave higher 
charmonium. If Z(3930) is a good candidate of Xc2&P) char- 
monium as claimed by Belle and BaBar, we must explain why 
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FIG. 2: (color online). Our best fit (red histogram) to the cos#* 
distribution of 77 — > DD. Here, we adopt the same parameters as 
those used in Fig. 1. Blue dots with error bars are the Belle's data 
[5], The cyan histogram is the background measured by Belle [5], 
while the grey solid curve is our best fit to this background. 



the signal of^ c o(2P) is missing in the measured DD invariant 
mass spectrum of yy — > DD [5, 6] because the theoretical 
study of P-wave charmonia indicated that the mass of Xc-o(2P) 
is close to that of Xci(2P) [8, 10] and Xco(2P) dominantly de- 
cays into DD [9]. In this letter, we have proposed that the 
observed structure Z(3930) should be formed by two P-wave 
higher charmonia ^ e o(2P) and xaQP). To examine whether 
this conjecture is reasonable, we have performed the analysis 
of the DD invariant mass spectrum and cos 0* distribution of 
yy — > DD. Our study has illustrated that the experimental data 
of these quantities can be well reproduced when considering 
both Xco{2P) and Xci{2P) in the yy — > DD process, which 
supports our conjecture of Z(3930) structure. Our fit to the 
Belle's data [5] gives the resonance parameters as, 

Mass(GeV) Width (Me V) 



Xco(2P) : 
Xa(2P) : 



3.92 
3.94 



8.07 
11.98 



which provide the valuable information on Xco(2P) and 
Xdi^-P). In Ref. [7], Belle announced the observation of 
a narrow enhancement X(3915) with mass 3915 ± 3(stat) ± 
2(syst) MeV and width 17 + 10(stat) ± 3(syst) MeV in yy -> 
J/if/u). Very recently, BaBar also confirmed the observation of 
X(3915) [12], where its mass and width are 3919.4+2.2(stat)+ 
1 .6(syst) Me V and 13+6(stat) + 3(syst) MeV, respectively. We 
notice that the measured mass and width of X(3915) are in 
good agreement with the extracted resonance parameters cor- 
responding to^ e o(2f), where X(3915) was once explained as 
the candidate ofxco(2P) [9]. As indicated in Ref. [12], BaBar 
is carrying out the spin-parity analysis of X(3915), which will 
be helpful to clarify whether it is suitable to explain X(3915) 
as f-wave higher charmonium^ c o(2P). 



It is obvious that our observations presented here would be 
the smoking gun for the study of f-wave higher charmonia. 
Beside testing their resonance parameters by different phe- 
nomenological models, experimental examination of whether 
there exist Xco(2P) and Xci(2P) in the yy — > DD data will be 
an important and interesting research topic in future experi- 
ments, especially in Belle, BaBar, forthcoming Bellell [13] 
and SuperB [14]. We expect further theoretical and experi- 
mental studies on these two P-wave higher charmonia. 

Note added in proof : After completing our paper, there ap- 
pears a paper (arXiv: 1207.2651), in which the BaBar collab- 
oration claims that they have confirmed the existence of the 
charmonium-like resonance X(3915) decaying to J/i//a> with 
a spin-parity assignment J p = + [15], i.e., the identification 
of the signal as due to the^ f o(2P) that we have claimed in this 
letter as well as in Ref. [9]. 
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